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SUMMARY 


The basic model of Lindhard and Scharff, known as the local plasma model, is 
utilized to study the effects on stopping power of the chemical and physical state of 
the medium. Unlike previous work with the local plasma model, in which individual 
electron shifts in the plasma frequency were estimated empirically, the Pines correc- 
tion derived for a degenerate Fermi gas is shown herein to provide a reasonable esti- 
mate, even on the atomic scale. Thus, the model is moved to a complete theoretical 
base requiring no empirical adjustments, as characteristic of past applications. The 
principal remaining error is in the overestimation of the low-energy absorption prop- 
erties that are characteristic of the plasma model in the region of the atomic dis- 
crete spectrum, although higher-energy phenomena are accurately represented, and even 
excitation-to-ionization ratios are given to fair accuracy. Mean excitation energies 
for covalent-bonded gases and solids, for ionic gases and crystals, and for metals 
are calculated using first-order models of the bonded states, for which reasonable 
agreement with the recently evaluated data of Seltzer and Berger is obtained. Hence, 
the methods described herein allow reasonable estimates of mean excitation energy for 
any physical-chemical combination of material media for stopping-power applications. 


INTRODUCTION 

The impact of energetic ions with atoms and molecules is a fundamental physical 
phenomenon with numerous applications in areas such as astrophysics, medical physics, 
high-power lasers, and plasma fusion devices, and is of special interest for radia- 
tion protection of people and materials in space. Most of what is known in modern 
physics is ultimately related to observations of charged-particle impact phenomena, 
by which particle character is studied, starting with the early works of Rutherford 
using radioactive sources. To advance our knowledge of such processes is both timely 
and technologically important. The historical context in which present theoretical 
developments are to be understood is discussed in the following section. 


HISTORICAL PERSPECTIVE 

It was recognized early in the classical treatment of charged particle slowing 
down that the f ree-electron, long-range Coulomb interaction leads to divergencies in 
the energy-loss rate. These divergencies indicate that there is a need for a long- 
range saturation effect. The saturation in gases was discussed by Bohr (ref. 1) in 
terms of Ehrenfest's principle. Bohr proposed that the saturation in gases is caused 
by the bonding of the electrons. To effect energy transfer, the interaction time 
T = b/v (where b is the impact parameter and v the ion velocity) must be short 
coirpared with the oscillating period of the bonded electron. (A list of symbols and 
abbreviations used in this paper appears after the references.) Hence, the adiabatic 
long-range collisions provide the necessary saturation, and an upper limit is estab- 
lished for the effective impact parameters. Most of our modern understanding stems 



from Bethe's detailed quantum theory (ref. 2) based on the Born approximation 
Stopping power for gaseous media with this approximation is given by 
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where Z^ is the projectile charge, n is the number of targets per unit volume, 

Z^ is the number of electrons per target, m is the electron mass, v is the pro- 
jectile velocity, 3 = v/c, c is the velocity of light, C is the velocity- 
dependent shell-correction term (ref. 3), and I 2 is the mean excitation energy 
given by solving 
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where f^ is the electric dipole oscillator strength of the target and is the 

corresponding excitation energy. The sum in equation (2) includes discrete and 
continuum levels. Empirically, it was observed that molecular stopping power is 
reasonably approximated by the sum of the corresponding empirically derived ’’atomic” 
stopping powers (ref. 4). Equations (1) and (2) imply 
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where Z and I pertain to the molecule, Zj and are the corresponding atomic 

values, and represents the stoichiometric coefficients. This additivity rule, 

given by equation (3), is called Bragg’s rule. 

Sources of deviations from Bragg’s additivity rule for molecules and the con- 
densed phase are discussed by Platzman (ref. 5). Aside from shifts in excitation 
energies and adjustments in line strengths as a result of molecular bonding, new 
terms in the stopping power are caused by the coupling of vibrational and rotational 
modes. Additionally, in the condensed phase, some discrete transitions are moved 
into the continuum, and collective modes among valence electrons in adjacent atoms 
produce new terms to be dealt with in the absorption spectrum. Platzman proposed 
that the experimentally observed additivity rule may not show that molecular stopping 
power is the sum of atomic processes but rather that it demonstrates that molecular 
bond shifts for covalent-bonded molecules are relatively independent of the molecular 
combination. On the basis of such arguments, Platzman suggested ionic-bonded sub- 
stances as the ones which should be studied as a rigid test of the additivity rule 
because of the radical difference in bonding type. He further estimated that ionic- 
bond shifts could change the stopping power by as much as 50 percent. 

Among the early indicators of the violation of the Bragg rule was the calcula- 
tion of 15 eV for the mean excitation energy of atomic hydrogen (using eq. (2) with 
the exactly known oscillator strengths and excitation levels) compared with a rather 
firmly established experimental value for molecular hydrogen of about 18 eV. Since 
accurate values of atomic mean excitation energies have been calculated for numerous 
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elements by Inokuti and co-workers (refs. 6 to 8) for the purpose of evaluating chem- 
ical bonding effects in molecules, empirical values have been substantially perturbed 
by effects of the chemical bonds. Although the mean excitation energy for gas mole- 
cules could be evaluated in principle from equation (2), the lack of knowledge of the 
excitation levels and corresponding oscillator strengths is the main hindrance. 

It was suggested by Dalgarno (ref. 9) that the oscillator strength distributions 
could be determined empirically from the photoabsorption spectra (aside from experi- 
mental uncertainty). Much of these data are obtained by energy-loss experiments by 
electron impact scattering at forward angles. Values of mean excitation energy for a 
number of simple molecules have in this way been estimated, and have demonstrated the 
shift in atomic values caused by chemical bonding (refs. 10 and 11). 

Theoretical calculation of mean excitation energies is hindered by the diffi- 
culty of solving for the complete excitation spectrum of complex quantum systems. 
Dalgarno was able to simplify the calculation by introducing a generalized function, 
which is related to the excitation spectrum as follows: 
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However, this function can be evaluated without explicitly forming the indicated 
sum. Thus, Dalgarno (ref. 12) was able to reduce equation (4) to 
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where i(; is the ground-state wave function, is the corresponding energy, o) 

is an energy eigenvalue, and X is the corresponding eigenvector. Chan and Dalgarno 
(ref. 13) calculated I as 42 eV for helium and Kamikawai et al. (ref. 14) calcu- 
lated 18.2 eV for molecular hydrogen by the same method. These values are in excel- 
lent agreement with experiments. 

Simultaneous with the development of the microscopic theory of stopping power 
was the macroscopic electrodynamic description of energy loss as required for the 
description of the long-range part of the interaction in condensed phase. This is 
because the interaction is simultaneous among many constituents. The slowing down is 
through the force exerted on the passing particle by the electric field induced in 
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the medium by the passage (ref* 15). It is customary to assume that the electric 
displacement vector is linearly related to the time-varying electric field as 
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for which the dielectric constant is 


e(a>) = 1 + / g{x) exp(icoT) dr (8) 

0 


The short-range collisions are still treated by Bethe theory with the result for 
total stopping power (see ref. 16 for details) of 
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where 6 is a density-effect correction applicable at high energies 
> 1/e (0)) . Also, 
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where Im(Z) denotes the imaginary part of Z. A result of dispersion theory is 
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where f(o)) is the dipole oscillator strength per unit cell of the medium, so that 

oo 
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which reduces to the usual Bethe expression (eq. (2)) in a sparse gas for which 
e (co) « 1 . 

If the long-range saturation effect is in terms of adiabatic limits for a gas 
and in terms of the medium polarization response for condensed dielectrics, the 
saturation effect for a f ree-electron gas is related to the tendency of a neutral 
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plasma to screen a local charge imbalance at large distances (ref, 17), The dielec- 
tric function of a f ree-electron gas is derived by Lindhard (ref. 18) and applied to 
the stopping-power problem for a classical electron gas and for the interaction-free 
Sommerfeld electron gas model. For a f ree-electron gas at rest, Lindhard arrives at 
the equation 
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where p 


is the electron density and w is the classical plasma frequency given by 


0) 


2 


P 


2 

4ire^ 

m 


P 


(14) 


Strictly speaking, equation (13) applies only when the electron gas is at rest, but 
it also applies in the limit of high projectile velocity compared with the average 
motion of the electrons. 

A discovery which paralleled the Lindhard investigations was made by Bohm and 
Pines (refs. 19 to 21), in which collective long-range interactions in a quantum 
electron gas were separated from individual electron motion through a canonical 
transformation, after which the normal coordinates of collective oscillation appear. 
This separation of the Hamiltonian into collective and individual electron motion is 
accomplished because of the effective screening of the Coulomb fields of individual 
electrons for distances greater than the screening distance For collective 

motion to give a major contribution to the Hamiltonian, the individual electron wave- 
length must be greater than X • Bohm and Pines (ref. 21) found the average collec- 
• c 

tive plasma frequency to be 
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where X is the average electron separation and x is the ratio of the average 
electron^wave length to the screening distance. Pines (ref. 22) suggests that the 
screening parameter x should be chosen to minimize the electron long-range correla- 
tion energy (that is, the electronic Coulomb energy), which, for plane-wave states 
appropriate to their degenerate electron gas, is given by 
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Pines (ref. 22) derived the stopping power in this degenerate electron gas and showed 
that the usual classical plasma frequency w is replaced by <( 0 >, which includes 
corrections for individual electron motion. ^ 
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A rather bold suggestion was made by Lindhard and Scharff (ref. 23) that equa- 
tion (13) could be applied on the atomic scale if the appropriate average over the 
atomic electron density was made. They further suggested that the effects of indi- 
vidual bonding of the electrons in their atomic orbitals could be incorporated 
through the added factor y « f2 as 
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From equation (17), the mean excitation enerqy is given by 
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Lindhard and Scharff estimated the mean excitation energy for atomic Hg as 768 eV 
compared with »800 eV from experiments. For He, they got 37 eV compared with 3 5 eV 
from quoted ejq>eriments. Ihey further approximated molecular hydrogen by taking the 
effective charge to be Z = 1 .2 and obtained 16 eV. 

Following this initial success of treating atoms as localized electron plasmas, 
Lindhard and Winther (ref. 24) extended equation (17) by using the more general 
velocity-dependent dielectric function derived by Lindhard (ref. 18), and demon- 
strated the ability of the Lindhard theory to predict tight bonding corrections of 
similar character to those of Walske (ref. 25) in connection with the Bethe theory. 

Chu and Powers (ref. 26) made extensive use of the work of Lindhard and Scharff 
(ref. 23) to demonstrate oscillations in the mean excitation energy. This work 

gave rise to corresponding oscillations in stopping power from which periodic 

variations are associated with the atomic shell structure (ref. 26). The more 
detailed calculations of Rousseau et al. (ref. 27) utilized the velocity-dependent 
Lindhard-Winther theory (ref. 24) and Bonderup*s simplified form of the Lindhard 
theory (ref. 28) and show good agreement with 2-MeV alpha -particle stopping-power 
data (ref. 29). Throughout these efforts, the parameter y is taken as the square 
root of 2, as suggested by Lindhard and Scharff (ref. 23). 

Chu et al. (ref. 30), using the theory of Lindhard and Winther in which individ- 
ual electron corrections to the local collective excitation were treated empirically 
by taking Y as an adjustable parameter, evaluated the aggregation effects for con- 
densed noble gases and metals. The condensed-gas calculations determined electron 
densities according to atomic Hartree-Fock densities, including overlap from the 
nearest neighbors in the condensed phase. Metallic wave functions were taken from 
the muffin-tin calculations of Moruzzi et al. (ref. 31). In most cases, y was in 
the range from 1.2 to 1.3. (See ref. 32.) 

As noted by Dehmer et al. (ref. 6) equation (18) may be rewritten as 
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from which can be obtained 


f (to) 



p(r) 


( 20 ) 


where 6(x) is the delta function. It is seen from equation (20) that, in the local 
plasma approximation, the volume of plasma with cutoff frequency y^p ~ ^ approxi- 
mates the total oscillator strength of the system at frequency (o. No exact equiva- 
lence is implied between the oscillator frequency distribution given by equation (20) 
and the oscillator frequency distribution of a quantum system. (Tlriis is true because 
eq. (20) exhibits a continuous spectrum, although quantvun systems generally exhibit a 
series of poles associated with the discrete quantum levels as well as a continuum at 
higher frequencies. ) Some insight may be gained by conparing dispersion relations 
for atomic systems with those for a related plasma. Ihe dispersion relation for a 
classical plasma is given by the dielectric constant e:(a)) as 
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where (jo^ is the usual plasma frequency and equation (21) results from the plasma 
conductivity (ref. 33). Indeed, the same pole appears in metals as the result of the 
conduction electrons from which the characteristic optical properties of metals 
derive (refs. 33 and 34). The more general dispersion relation, derived from equa- 
tions (8) and (11), is 
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where P denotes the principal value at the singularity. In atomic systems, the 

oscillator strengths are broadly separated in frequencies according to shells; the 

outer shells appear at the lowest frequencies, and the innermost shell appears at the 

highest frequencies. The lack of oscillator strength at frequencies between shells 

results in large gaps in the spectrum. Let w be a frequency in the broad gap 

between two successive shells - the first centered at o) and the second at w . 

1 2 

Then the dispersion relation (eq. (22)) becomes 
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where 
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so that is the plasma frequency associated with the electrons of the outermost 

shell. Although equations (23) and (24) provide motivation for using the local 
plasma approximation (eq. (20) ), there is plenty of room for a more complete under- 
standing as to why the model works as well as it does in practical calculations 
(ref. 35). 

In previous investigations, we considered the use of the local plasma model to 
evaluate molecular bonding effects on the mean excitation energy of molecules of 
covalently bonded hydrogen and carbon (ref. 35) as well as ionic crystals and gases 
(ref. 36), in which quite sensible corrections to the usual Bragg rule were obtained. 
The chemical bond shifts were unambiguously defined in terms of atomic integrals and 
molecular parameters* In the usual implementation of the local plasma model 
(eq. (18)), y corrects for a shift in the local plasma frequency caused by individ- 
ual electron effects. Lindhard and Scharff (ref. 23) suggest y = f2; however, 
y » 1 .2 yields atomic mean excitation energies from the local plasma model in better 
agreement with the accurate atomic values calculated by Dehmer et al. (ref. 6). The 
fact that the larger value (y = ^) gives better agreement with empirical data sug- 
gests that this larger value corrects (in addition to individual electron shifts) for 
the chemical shifts as well. Such chemical shifts were estimated separately for 
covalent and ionic bonds in references 35 and 36. 

Encouraged by the smallness (<30 percent) of the empirical individual electron 
corrections to the collective plasma frequency (refs. 32, 35, and 36), a calculation 
(ref. 37) in which individual electron shifts were estimated according to the theory 
for plane-wave states in an extended plasma, as calculated by Pines (ref. 22), yields 
results that are in good agreement with Dehmer et al. (ref. 6). Consequently, the 
local plasma model is placed on a parameter-free basis in which chemical shifts are 
determined from atomic molecular parameters alone, and effects of individual electron 
motion are evaluated in terms of the Pines correction, the combined effects of which 
are on the order of the plasma frequency shift of y ” suggested by Lindhard and 

Scharff. 

The Pines correction makes a remarkable improvement in the prediction of the 
local plasma model, and further adjustments in the theory, to account for the plasma 
frequency shifts resulting from the atomic shell structure, should bring the model 
into accurate predictive capability. To further elucidate the relationship between 
the local plasma model and the more exact quantum treatment of bonded systems, 
related quantities of both theories in the case of one- and two-electron systems are 
examined in the following section of this report. Atomic mean excitation energies 
and straggling parameters, based on the local plasma model, are compared with accu- 
rate calculations of Inokuti et al, (refs. 6 to 8) in the section entitled "Stopping 
and Straggling Parameters of Atoms." The use of the Gordon-Kim electron gas model of 
molecular bonding (ref, 38) to determine the effects of covalent chemical bond shifts 
on mean excitation energy for elements of the first two rows is presented in the 
section entitled "Covalent-Bond Effects." Calculations of mean excitation energies 
of ionic-bonded substances are discussed in the section entitled "Ionic-Bond 
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Effects," and the mean excitation energies of metals are discussed in the section 
entitled "Metal lie -Bond Effects." 


EXCITATION SPECTRA OF ONE- AND TWO-ELECTRON SYSTEMS 

The hydrogen atomic excitation spectrum in the dipole approximation is well 
known as 
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where n is the principal quantum number, R is Rydberg's constant, oj is given 
by 


boj = R(1 - 1/n ) 
n ' 


(26) 


and 


Rk = fi03 - R 


(27) 


The corresponding spectrum for the local plasma model (eq« (20)) is given as 
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Where fito^ = 55.12 eV. The cumulative oscillator strength 


(0 

F(o)) = / f(aj') dW (29) 
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is shown in figure 1 for each of the two models. Similarly, the excitation spectrum 
of the helium atom has been evaluated for screened wave functions and is shown in 
figure 1 . The fractional excitations of the two models never differ by more than 
“15 percent above the excitation threshold. As noted by Dehmer et al. (ref. 6), the 
main error in the local plasma model is the contribution to absorption below excita- 
tion threshold all the way down to zero* This error is also evident in the energy 
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moments of the plasma model. The moments of the energy spectrum for the hydrogen 
atom are shown in figure 2 , where 


<(E/R)"'> = 



f(o)) do) 


( 30 ) 


and m is a continuous parameter. The low-frequency contributions associated with 
the local plasma model cause a divergence in equation (30) at m « -2 which is not 
present in the quantum system. Atomic polarizability and the low-frequency refrac- 
tive index are affected the most. Other atomic properties, such as the total inelas- 
tic cross section, the mean excitation energy, the straggling parameter, and the mean 
electronic kinetic energy, are reasonably represented by the plasma model. Also 
shown in figure 2 are results, including the Pines correction (ref. 22) to the plasma 
frequency, which indicate substantial improvement in the prediction of atomic proper- 
ties, although low-energy atomic properties are still beyond the scope of the model. 

The plasma model is expected to be more accurate as more electrons are added to 
the system. This occurs in two ways, as seen in figure 1. First, a greater contri- 
bution comes from the continuum, which is most like the plasma. Second, the excita- 
tion thresholds shift to relatively lower energies and fill in the low-frequency 
region, for which the plasma model normally tends to error. A considerable improve- 
ment in the energy moments of helium for the local plasma approximation are clearly 
shown in figure 3. 

The moments of the excitation spectrum of have been evaluated empirically 

by using experimental oscillator strengths (ref. 39) and theoretically (ref. 14) 
using the Dalgarno sum rules of reference 12 (eqs. (4) to (6)). These are compared 
in figure 4 with an "atomic” approximation to taken as a generalization of 

Bragg's rule (ref. 11). Also shown in figure 4 are values for calculated using 

the local plasma model with the Pines correction and with the Gordon-Kim model of the 
molecular wave functions (ref. 38) given as 


pR (r) = pR(r) + pR(r - R) (31) 

2 

where atomic hydrogen electron density and R is the displacement 

vector of length between the two centers. It is clear from figure 4 that, 

even with the simple Gordon-Kim approximation, the plasma model is a considerable 
improvement over the Bragg rule, except for the lowest-energy molecular properties 
(i.e., m < -0.5). Figure 4 also shows that the Gordon-Kim approximation introduces 
minor errors compared with the inherent limitations of the local plasma model. 

The mean excitation energy for stopping power may likewise be evaluated. Atomic 
and molecular hydrogen are compared in table I with a recent compilation of experi- 
mental data (ref. 40). Quite reasonable estimates of atomic and molecular properties 
of importance to ionizing radiation are obtained by this local plasma model if the 
Pines correction is included. Optical and other low-frequency properties, however, 
are poorly represented. The plasma model should become more accurate for more com- 
plex many-electron systems, especially those in which the optical properties are more 
in line with those predicted by the plasma model. 

With the present results, it is now clear what approach should be taken to 
improve the plasma model applications. Clearly, a correction factor similar to that 
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of Pines should be introduced to suppress absorption below excitation threshold and, 
correspondingly, to enhance frequencies just above threshold. A number of possibili- 
ties are open to implement such a correction, which would appear as a first-order 
quantum correction for the discrete spectrum. Preliminary work by Walecka (ref. 41) 
on the study of collective atomic oscillations may be a starting point for further 
development. 


STOPPING AND STRAGGLING PARAMETERS OF ATOMS 

In this section, parameters are considered for atoms associated with the stop- 
ping of charged particles and fluctuations in their energy transfer. The energy 
moment is 


S(m) = <(E/R)®> 


(32) 


and the related quantity is 


L(m) 


dS(m) 

dm 


(33) 


In terms of these quantities, the mean excitation energy is 


In I = L(0)/S(0) 


(34) 


and the straggling parameter related to fluctuations in energy loss is 


In A = L(1 )/S(l ) 


(35) 


The mean excitation energy (eq. (34)) has been evaluated in the context of the local 
plasma model and is compared in figure 5 with the values computed by Inokuti and 
co-workers (refs. 6 to 8) for atoms through krypton. Hartree-Fock wave functions 
(ref. 42) have been used through neon and sodium through krypton are represented by 
screened wave functions (ref. 43). 

The values for the straggling parameter were similarly evaluated and are com- 
pared in figure 6 with the values obtained by Inokuti and co-workers. Also shown are 
values for noble gases coitpiled by Inokuti et al. (refs. 6 to 8) and values obtained 
by Zeiss et al. (ref. 11). The present values tend to be about 25 percent low at 
Z « 36, with improvements at lower values of Z, which may be caused by the lack of 
shell structure corrections in the plasma frequencies of the K and L shells. 

It is clear from these atomic calculations that the plasma model with the Pines 
correction generally provides good results for mean excitation energy and reasonable 
estimates for the straggling parameter. Although Hartree-Fock wave functions are 
required for low atomic numbers, reasonable results are obtained using screened wave 
functions for atoms heavier than argon. It is mainly the low-energy atomic 


properties which require improvements beyond the Pines correction. *Ihese properties 
emphasize the need for a first-order quantum correction to the atomic structure. 


COVALENT-BOND EFFECTS 

Early experimental work with ionization energy loss was conducted in covalent- 
bonded gases (also noble gases), from which Bragg's rule was derived. Although more 
recent experimental work, beginning with Thompson (ref. 44), has shown systematic 
variation from Bragg's rule, such rules still seem appropriate for fixed molecular 
structures (refs. 45 and 46). As the result of the theoretical efforts of Inokuti 
and co-workers (refs. 6 to 8) , it is clear that chemical-bond shifts in the mean 
excitation energy have occurred, and, as suggested by Platzman (ref. 5), all covalent 
shifts are of similar magnitude. 

In any molecular dynamic calculation, there is a trade-off between model accu- 
racy and computational efficiency. As pertains to the radiolysis of large molecular 
structures, the most useful model is the lowest order possible. It is clear that the 
use of self-consistent field methods to determine molecular wave functions would 
seriously limit the ability to study systems of practical interest. Considering the 
relative success of the Gordon-Kim electron gas model of molecular bonding (refs. 38, 
47, and 48) , a simple method for the calculation of chemical-bond effects on the mean 
excitation energies is suggested. As suggested by Gordon and Kim, the molecular 
electron density as a superposition of the unperturbed atomic states is given by 


p(r ) = (r ) + 2^ 


(36) 


for diatomic molecules. There is an obvious generalization of equation (36) for the 
polyatomic case. Whereas Gordon and Kim used equation (36) to calculate the molecu- 
lar potential (see refs. 47 and 48 for ionic and covalent applications) from which 
2 is theoretically obtained, is taken here from observed experimental bond 

distances. Substituting equation (36) into equations (14) and (18) and reducing 
results in 
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where I^ and I 2 are the corresponding atomic values, which are accurately known 
(refs. 6 to 8), The chemical bonding correction is generally 
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Generalizing equation (37) yields 
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where the sum over j includes every bond in which Z is attached in the molecule. 
Correction factors have been calculated (ref, 35) for hydrogen and carbon molecules 
with the bond parameters in table II, Carbon sp^ hybrid orbital wave functions were 
used in these calculations, although values were only slightly different. The 

tetrahedral orbitals were spherically symmetrical in their electron densities. 
Therefore, spherical symmetry was assumed throughout subsequent calculations. 

Recommended values of mean excitation energies (ref , 40 ) are compared in 
table III with theoretical values calculated using atomic mean excitation energies 
from Dehmer et al, (ref. 6) with the bond corrections in table II, Bragg’s rule is 

also used with the atomic values of Dehmer et al, for comparison. Although the theo- 

retical values are within 4 percent of the experimental and empirical values, the 
Bragg’s rule values are from 17 to 21 percent low, indicating a substantial adjust- 
ment as the result of chemical bonding. 

Mean excitation energies have been calculated for covalent gases of the first 
two rows using the local plasma model and the Pines correction. Results are compared 

with empirical values (ref, 40) in table IV, Corresponding values for covalent sol- 

ids are shown in table V. 

Moments for the ^2 molecule using the plasma model are compared in figure 7 with 
values calculated from the oscillator strengths compiled by Dalgarno et al. 

(ref, 49), As can be seen, good agreement between the present simple plasma model 
calculations and the oscillator strength distribution of reference 49 is obtained 
except for the lowest-f requency phenomena. 


IONIC-BOND EFFECTS 

Although covalent bond shifts were found to be relatively small corrections to 
atomic values, such a separation as in equation (37) in terms of neutral atomic val- 
ues are not possible for ionic bonds. Using the Gordon-Kim model electron density of 
the partial ionic (diatomic) system, 

A rJ 

where and B^“^^ refer to partially ionic states of the two constituents, 

^AB their nuclear separation, and p is the partial ionic fraction. The elec- 
tron density of a partial ionic atom in equation (40) is 




where Pj.(r) is the electron density of the neutral atom and p (r) is the elec- 

A^ 

tron density of the atomic ion. With the aid of equations (40) and (41), shifts in 
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the mean excitation energy caused by ionic and covalent effects can be evaluated. As 
shown in reference 36, 


Z ln{I) = 


with 
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(43) 


where y is the Pines correction given by equation (15) or estimated empirically as 
in reference 46. Mean excitation energies for various stages of ionization calcu* 
lated using the Pines correction and the atomic wave functions of Clementi and Roetti 
(ref. 42) are shown in figure 8. In addition to the the ionic-bond shifts, there are 
shifts caused by covalent-like character, as given by 


In 1 
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Mean excitation energies for partial ionic-bonded substances are shown in table VI 
with the corresponding bond parameters used in the model. Also shown are values for 
a pure covalent bond and Bragg values using the neutral atomic mean excitation ener- 
gies of Dehmer et al. (ref. 6) as well as Bragg values of the corresponding partial 
ionic states. The ionic-bond fractions are taken from Pauling (ref. 50) as experi- 
mental data for HF and LiH. Bond lengths are for ionic crystals except for the HF 
gas. Atomic mean excitation energies are shown for partial ionic states in figure 8, 
and differ from values in reference 36 because of the Pines correction. 

It is clear from table VI that the main contribution to corrections to the Bragg 
rule is the adjustment from atomic neutral to atomic ion mean excitation energies as 
proposed by Platzman (ref. 5). Indeed, when there is little difference between the 
usual Bragg value and the partially ionic Bragg value, the covalent value is in near 
agreement with the predicted value of I for HF and LiH in the table. For LiF, the 
relatively large adjustment from the usual Bragg value (81.6) to the partially ionic 
Bragg value (92.6) leaves a large difference between the covalent value (83.4) and 
the predicted value of I (93.6). The adjustment of the ionic-bond shift caused by 
the covalent-like character for LiF is 1 eV compared with adjustments of the neutral 
states caused by the pure covalent bond of 1.8 eV. This comparison shows the greater 
role of the Coulomb attraction in forming the bond of the ionic molecules relative to 
the two-electron interaction in forming the covalent bond. 

Calculated mean excitation energies for ionic crystals using the Pines correc- 
tion are shown in table VII, along with recommended values (ref, 40). The crystal 
parameter and fractional ionic charge have been taken from Pauling (ref. 50). The 
LiF value is the only one with an experimental basis (ref, 36), 



METALLIC-BOND EFFECTS 


The first approach to metals in this paper is similar to that taken by Chu et 
al. (ref* 30) ^ in which they employed the muffin-tin wave functions (ref. 31) and 
stopping-power theory according to Lindhard and Winther (ref. 24)* Individual elec- 
tron corrections to the local plasma frequency are treated empirically through an 
adjustable parameter y (See table I of (ref. 32) and related discussion.) Unlike 
this previous work, the present work includes estimates of shifts in the plasma fre- 
quency according to the Pines correction in equation (15) and is in that sense com- 
pletely deterministic. 

The metallic wave functions for lithium metal approximated by the Wigner-Seitz 
model (ref. 51) are considered first. In deriving these wave functions, the lithium 
ion core potential was taken from the screened wave functions of Clementi and 
Raimondi (ref. 43), and the calculated crystal -valence wave functions (aside from 
normalization) were found to be a slight perturbation (mainly due to boundary condi- 
tions) of the free hydrogenic (2s) orbital inside the Wigner-Seitz sphere. The final 
crystal wave functions used were constructed from unperturbed Hartree-Fock orbitals 
(ref. 42) in the core region with a small perturbation outside the core. This per- 
turbation matched the boundary conditions on the surface of the Wigner-Seitz sphere 
(ref. 51). This was followed by normalization of the valence-shell wave functions 
(to make the valence electron density add up to give the correct number of valence 
electrons). These wave functions are quite similar to the muffin-tin model (ref. 31) 
and yield mean excitation energies in substantial agreement with Ziegler (ref. 32) 
when Y is taken as his empirical value. The mean excitation energies for metals of 
the second and third row using Wigner-Seitz wave functions (treating all valence 
electrons as spatially equivalent) and the Pines correction are compared in table 
VIII with empirical values (ref. 40) • 

The present results clearly demonstrate that the effects of the metallic bond in 
lithium and beryllium are large and are mainly the result of collective oscillations 
in the f ree-electron gas formed by the valence electrons. Although similar good 
agreement should be expected for sodium and magnesium, it is emphasized here only 
that these empirical values are interpolations without an experimental basis, and 
smaller empirical values more in line with the present results should not be elimi- 
nated. The small value predicted for aluminum (149 eV) is in doubt, as the empirical 
value (166 ± 4 eV) is based on one of the most experimentally studied quantities 
since aluminum served as a standard in stopping-power experiments for many years. 

The fault could well lie in the use of the Wigner-Seitz model for group III metals. 

It is well-known that the success of the Wigner-Seitz theory rests mainly on applica- 
tion to alkali metals. Although some hope for application to group II metals exists, 
treating the three valence electrons of group III as spatially equivalent is clearly 
in error. Correction to metals from an alternate model, proposed by Pauling 
(ref. 50) for metallic orbitals and implemented here in simplified fashion, are 
considered next. 

In X-ray diffraction experiments, even beryllium metal shows a considerable 
degree of covalent quality, as suspected from bulk material properties (ref. 52). In 
this view, a model is considered in which the valence-bond effects can be included 
explicitly. In the spirit of the Pauling valence-bond theory (ref. 50) and the 
Gordon-Kim model (ref. 38) of valence bonding, the electron density about the ion 
cores is assumed to be a superposition of partial ionic core states among nearest 
core neighbors. Additional contributions from next -nearest neighbors are assumed to 
add to the electron continuum states in a manner analogous to the Pauling 
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unsynchronized resonances in lithium crystals (ref. 50), Ihe electron density of the 
partially ionic core charge p is 


P 

A 


(+p) 


(r) 



(?) 



(45) 


is the electron density of the 


where v is the number of valence electrons, 
atomic neutral state, and p^(+^)(3^) is the electron density of the valence-stripped 


ion cores. We have used the observation by Slater that radial wave functions of the 
L shell are nearly the same for both values of as a result of exchange interac- 

tion between the (2s) and (Is) orbitals. The same is true for the M shell. In the 
present calculation, each metal ion core has been placed into a Wigner-Seitz cell and 
the electron density from nearest neighbors has been approximated by reflecting the 
exterior core density function across the cell boundary. The continuum electron 
density is then taken as 


p = [p + (v - p)6](3/47i:r ^) 
e s 


(46) 


where 6 is the next -nearest-neighbor contribution to the continuum. The value of 
6 is determined by requiring a full cortplement of v valence electrons per cell. 

The resultant electron density p(r) was used to calculate the local plasma fre- 
quency and mean excitation energy per cell. The Pines correction was used for indi- 
vidual particle shifts. The radii r^ for the Wigner-Seitz cell are given in 
table VIII. The ion-core wave functions were calculated from the Hartree-Fock wave 
functions of dementi and Roetti (ref. 42). A slight dependence on the ion-core 
charge appears (ref. 53) in which there is some increase in mean excitation to 
I = 155 eV for aluminum. However, there are some unresolved questions, concerning 
periodicity at the cell boundaries, which leave the value of this model in doubt. 

The mean excitation energy for aluminum requires the reconsideration of the data 
on which it is based and the corresponding analysis. In an analysis by Andersen and 
Ziegler (ref. 54), 162 eV was assumed as the mean excitation energy for aluminum. A 
reduction in I to 1 50 eV results in a 3-percent increase in stopping power at 
1 MeV, which leaves it within the stated uncertainty limits of the Andersen and 
Ziegler parametric curves (ref. 54). These curves correspond to the uncertainty in 
the experimental data used in the analysis of reference 54. (See fig. 9.) Indeed, a 
number of authors have reported mean excitation energies for aluminum in line with 
the present results (refs. 55 to 60), although more recent analyses are higher. A 
recent study of aluminum optical properties indicates that a value of I several eV 
lower than 166 eV is not inconsistent with the empirical dielectric function 
(ref. 61). The shift of several eV is associated with polarization of the Al core 
by the valence electrons in their metallic orbitals. Such core polarization effects 
are not calculated in the present model. Furthermore, quantum corrections to K- and 
L-shell discrete spectrum may cause further small adjustments. In any case, the 
apparent discrepancy is due to the electronic wave functions used in the present 
calculation, to the inadequate treatment of corrections to the Bethe formula, from 
which I is extracted from the experimental data (see for example refs. 62 and 63), 
to quantum corrections, or to a combination of these. 

To further clarify the relationship between the mean excitation energy for alu- 
minum and experimental data, a band is shown in figure 9 which brackets the experi- 
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mental data in references 62 and 64 to 68 for proton energies between 0.5 and 
10 MeV. These energies are compared with the reduced stopping power calculated from 
the- Andersen and Ziegler empirical shell corrections (ref. 54). The older data of 
Kahn (ref. 64), which would have lowered the band considerably, were excluded from 
the figure. The mean excitation energies exhibited in the figure are the 167 eV used 
as input to Shiles et al. (ref. 61), the 162 eV determined by Andersen and Ziegler 
(ref. 54), the 155 eV estimated using one form of valence-bond theory, and the 149 eV 
calculated according to the present (simplified) Wigner-Seitz model. While it is not 
clear that the 167-eV curve is superior to the 149-eV curve, a modest shift in the 
empirical shell corrections can bring any of the four curves into an equally good fit 
to *the data. It is further emphasized that shell corrections are not exactly known 
nor in empirical analysis are shell corrections usually differentiated from other 
corrections to the Bethe formula (eq. (1)). 


DISCUSSION 

The present results are combined in figure 10 with the evaluated data of Seltzer 
and Berger (ref. 40). Care is taken where possible to model the same physical- 
chemical state. (See specific tables for details.) Results for free atoms (Hartree- 
Fock wave functions for Z < 1 0, and screened wave functions elsewhere) are also com- 
pared in figure 10 with the accurate atomic values of Dehmer et al. (ref. 6). It is 
clear that the trends in the first- and second-period elements are well approximated 
by the present application of the local plasma model, especially when the Pines cor- 
rection is applied. The present results are generally in fair agreement with the 
compilation and recommendations of Seltzer and Berger (ref. 40), although small dis- 
crepancies in the third period remain to be resolved. 

Perhaps the greatest criticism of the present application of the local plasma 
model calculations is the use of the Gordon-Kim approximation to the covalent-bonded 
wave functions. When the moments of the energy spectrum are considered, it is clear 
that the Gordon-Kim model approximately adjusts the excitation spectrum in the region 
of greatest importance to ionizing radiation, and appears no more in error than the 
basic plasma model in which it is used. (See fig. 4.) Of course, accurate use of 
the plasma model implies the necessary use of the Pines correction, as demonstrated 
for the hydrogen atom in figure 2 and used throughout the present calculations. 
Although the Pines correction produces marked improvements in the predictive capabil- 
ity of the model, further quantum corrections for the discrete spectrum would produce 
additional corrections and would hopefully remove most of the remaining error in the 
plasma model. Further improvement in electronic wave functions would be helpful in 
identifying the remaining corrections required for the plasma model. 


ACKNOWLEDGMENTS 

The authors wish to express their gratitude to Govind S. Khandelwal of Old 
Dominion University, Willard E. Meador of the Langley Research Center, Mito Inokuti 


17 


of the Argonne National Laboratory, and Stephen Seltzer of the National Bureau of 
Standards for helpful discussions and encouragement. Special thanks are offered to 
Stephen Seltzer, Martin Berger, and Mi to Inokuti for sharing unpublished data. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 3, 1 984 



REFERENCES 


1. Bohr^ N, ; On the Decrease of Velocity of Swiftly Moving Electrified Particles in 

Passing Through Matter. Philos. Mag. & J. Sci., ser. 6, vol. 30, no. 175, 

July 1915, pp. 581-612. 

2. Bethe, Von H.: Zur Theorie des Durchgangs schneller Korpuskularstrahlen durch 

Materie. Ann. Phys. , 5 Folge, Bd. 5, 1930, pp. 325-400. 

3. Walske, M. C.; and Bethe, H. A.: Asymptotic Formula for Stopping Power of 

K-Electrons. Phys. Rev., vol. 83, 1951, pp. 457-458. 

4. Bragg, W. H.; and Kleeman, R.: On the a Particles of Radium, and Iheir Loss of 

Range in Passing Through Various Atoms and Molecules. Philos. Mag. & J. Sci., 
ser. 6, vol. 10, no. 57, Sept. 1905, pp. 318-340. 

5. Platzman, Robert L. : Influences of Details of Electronic Binding on Penetration 

Phenomena, and the Penetration of Energetic Charged Particles Through Liquid 
Water. Symposium on Radiobiology - The Basic Aspects of Radiation Effects on 
Living Systems, James J. Nickson, ed. , John Wiley & Sons, Inc., c.195 2, 
pp. 139-176. 

6. Dehmer, J. L. ; Inokuti, Mitio; and Saxon, R. P. : Systematics of Moments of 

Dipole Oscillator-Strength Distributions for Atoms of the First and Second 
Row. Phys. Rev. A, third ser., vol. 12, no. 1, July 1975, pp. 102-121. 

7. Inokuti, Mitio; Baer, T.; and Dehmer, J. L. : Addendum: Systematics of Moments 

of Dipole Oscillator-Strength Distributions for Atoms in the First and Second 
Row. Phys. Rev. A, vol. 17, no. 3, Mar. 1978, pp. 1229-1231. 

8. Inokuti, M.; Dehmer, J. L.; Baer, T.; and Hanson, J. D. : Oscillator-Strength 

Moments, Stopping Powers, and Total Inelastic-Scattering Cross Sections of All 
Atoms Through Strontium. Phys. Rev. A, vol. 23, no. 1, Jan. 1981, pp. 95-109. 

9. Dalgarno, A.: The Stopping Powers of Atoms. Proc. Phys. Soc. (London), vol. 76, 

pt. 3, no. 489, Sept. 1, 1960, pp. 422-424. 

10. Zeiss, G. D.; and Meath, W. J. : The H^O-H^O Dispersion Energy Constant and the 

Dispersion of the Specific Refractivity of Dilute Water Vapour. Mol. Phys., 
vol. 30, no. 1, July 1975, pp. 161-169. 

11. Zeiss, G. D.; Meath, William J. ; MacDonald, J. C. F.; and Dawson, D. J. : Accu- 

rate Evaluation of Stopping and Straggling Mean Excitation Energies for N, O, 
^ 2 ' ^2' ^2' ^^ 3 ' ^2^' ^2^ Using Dipole Oscillator Strength Distribu- 

tions - A Test of the Validity of Bragg’s Rule. Radiat. Res., vol. 70, no. 2, 
May 1977, pp. 284-303. 

12. Dalgarno, A.: Sum Rules and Atomic Structure. Rev. Mod. Phys., vol. 35, no. 3, 

July 1963, pp. 522-524. 

13. Chan, Y. M.; and Dalgarno, A.: The Stopping Powers of Atoms and Molecules. 

Proc. Phys. Soc. (London), vol. 85, 1965, pp. 457-459. 


19 



14. Kamikawai, Ryotaro; Watanabe, Tsutoma; and Amemiya, Ayao: Mean Excitation Energy 

for Stopping Power of Molecular Hydrogen. Phys . Rev., second ser. , vol. 184, 
no. 2, Aug. 10, 1969, pp. 303-311 . 

15. Landau, L. D.; and Lifshitz, E. M.: Electrodynamics of Continuous Media. 

Addison-Wesley Pub. Co., Inc., 1960. 

16. Ahlen, Steven P. : Theoretical and Experimental Aspects of the Energy Loss of 

Relativistic Heavily Ionizing Particles. Rev. Mod. Phys., vol. 52, no. 1, 

Jan. 1980, pp. 121-173. 

17. Kramers, H. A.: The Stopping Power of a Metal for a-Particles. Physica, 

vol. 13, 1947, pp. 401-412. 

18. Lindhard, J. : On the Properties of a Gas of Charged Particles. Mat.-Fys. 

Medd. - K. Dan Vidensk. Selsk., vol. 28, no. 8, 1954. 

19. Bohm, David; and Pines, David: A Collective Description of Electron Interac- 

tions. I. Magnetic Interactions. Phys. Rev., second ser. , vol. 82, no. 5, 

June 1 , 1951, pp. 625-634. 

20. Pines, David; and Bohm, David: A Collective Description of Electron Interac- 

tions. II. Collective vs Individual Particle Aspects of the Interactions. 

Phys. Rev., second ser., vol. 85, no. 2, Jan. 15, 1952, pp. 338-353. 

21. Bohm, David; and Pines, David: A Collective Description of Electron Interac- 

tions. III. Coulomb Interactions in a Degenerate Electron Gas. Phys. Rev., 
second ser., vol. 92, no. 3, Nov. 1 , 1953, pp. 609-625. 

22. Pines, David: A Collective Description of Electron Interactions: IV. Electron 

Interaction in Metals. Phys. Rev., second ser., vol. 92, no. 3, Nov. 1 , 1953, 
pp. 626-636. 

2 3. Lindhard, Jens; and Scharff, Morten: Recent Developments in the Theory of Stop- 

ping Power. 1. Principles of the Statistical Method. Penetration of Charged 
Particles in Matter, Edwin Albrecht Uehling, ed. , Nucl. Sci. Ser. Rep. No. 29 
(Publ. 752), Natl. Acad. Sci. -Natl. Res. Counc. , 1960, pp. 49-55. 

24. Lindhard, J. ; and Winther, A.: Stopping Power of Electron Gas and Bquipartition 

Rule. Mat.-Fys. Medd. - K. Dan. Vidensk. Selsk., vol. 34, no. 4, 1964. 

25. Walske, M. C.: The Stopping Power of K-Electrons. Phys. Rev., second ser., 

vol. 88, no. 6, Dec. 1 5, 1952, pp. 1283-1289. 

26. Chu, W. K. ; and Powers, D. : Calculations of Mean Excitation Energy for All Ele- 

ments. Phys. Lett., vol. 40A, no. 1, June 19, 1972, pp. 23-24. 

27. Rousseau, C. C.; Chu, W. K. ; and Powers, D. : Calculations of Stopping Cross 

Sections for 0.8- to 2.0-MeV Alpha Particles. Phys. Rev. A, third ser., 
vol. 4, no. 3, Sept. 1971, pp. 1066-1070. 

28. Bonderup, E. : Stopping of Swift Protons Evaluated From Statistical Atomic Model. 

Mat.-Fys. Medd. - K. Dan. Vidensk. Selsk., vol. 35, no. 17, 1967. 


20 



29. Chu, W. K. ; and Powers^ D, : Alpha-Particle Stopping Cross Section in Solids From 

400 KeV to 2 MeV, Phys . Rev., second ser. , vol. 187, no. 2, Nov. 10, 1969, 
pp. 478-490. 

30. Chu, W. K. ; Moruzzi, V. L. ; and Ziegler, J. F.: Calculations of the Energy Loss 

of ^He Ions in Solid Elements. J. Appl. Phys., vol. 46, no. 7, July 1975, 
pp. 2817-2820. 

31. Moruzzi, V. L.; Janak, J. F.; and Williams, A. R.: Calculated Electronic Proper- 

ties of Metals. Pergamon Press, Inc., c.1978. 

32. Ziegler, J. F.: Ihe Stopping of Energetic Ions in Solids. Nucl. Instrum. & 

Methods, vol. 168, nos. 1-3, Jan. 1-15, 1980, pp. 17-24. 

33. Hubbard, J. : Ihe Dielectric Theory of Electronic Interactions in Solids. Proc. 

Phys. Soc. (London), vol. 68, pt. 11, no. 431A, Nov. 1, 1955, pp. 976-986. 

34. Frohlich, H.; and Pelzer, H.: Plasma Oscillations and Energy Loss of Charged 

Particles in Solids. Proc. Phys. Soc. (London), vol. 68, pt. 6, no. 426A, 

June 1, 1955, pp. 525-529. 

35. Wilson, J. W.; and Kamaratos, E.: Mean Excitation Energy for Molecules of Hydro- 

gen and Carbon. Phys. Lett., vol. 85A, no. 1, Sept. 7, 1981, pp. 27-29. 

36. Wilson, J. W.; Chang, C. K. ; Xu, Y. J. ; and Kamaratos, E.: Ionic Bond Effects on 

the Mean Excitation Energy for Stopping Power. J. Appl. Phys., vol. 53, no. 2, 

Feb. 1982, pp. 828-830. 

37. Wilson, J. W.; and Xu, Y. J. : Metallic Bond Effects on Mean Excitation Energies 

for Stopping Power. Phys. Lett., vol. 90A, no. 5, July 12, 1982, pp. 253-255. 

38. Gordon, Roy G. ; and Kim, Yung Sik: Theory for the Forces Between Closed-Shell 

Atoms and Molecules. J. Chem. Phys., vol. 56, no. 6, Mar. 15, 1972, 
pp. 3122-3133. 

39. Dalgarno, A.; and Williams, D. A.: Properties of the Hydrogen Molecule. Proc. 

Phys. Soc. (London), vol. 85, pt. 4, no. 546, Apr. 1965, pp. 685-689. 

40. Seltzer, S. M.; and Berger, M. J. : Evaluation of the Collision Stopping Power of 

Elements and Compounds for Electrons and Positrons. Int. J. Appl. Radiat. & 
Isot. , vol. 33, no. 11, Nov. 1982, pp. 1189-1218. 

41. Walecka, J. D. : Collective Excitations in Atoms. Phys. Lett., vol. 58A, no. 2, 

Aug. 9, 1 976, pp. 83-86. 

42. dementi, Ehrico? and Roetti, Carla: Roothaan-Hartree-Fock Atomic Wavef unctions. 

At. Data & Nucl. Data Tables, vol. 14, nos. 3-4, Sept. -Oct. 1974, pp. 177-478. 

43. dementi, E.; and Raimondi, D. L.: Atomic Screening Constants From SCF Func- 

tions. J. Chem. Phys., vol. 38, no. 11, June 1, 1963, pp. 2686-2689. 

44. Thonpson, Theos Jardin: Effect of Chemical Structure on Stopping Powers for 

High-Energy Protons. UCRL-1910, Univ. of California, Aug. 1 1 , 1952. 


21 


45. Lodhi, A. S.; and Powers/ D. : Energy Loss of a Particles in Gaseous C-H and 

C-H-F Compounds. Phys. Rev. A, vol. 10/ no. 6, Dec. 1974/ pp. 2131-2140. 

46. Neuwirth/ Wolfgang; Pietsch/ Werner; and Kreutz / Ronald: Chemical Influences on 

the Stopping Power. Nucl. Instrum. & Methods/ vol. 149, nos. 1-3/ Feb. 15- 
Mar. 1/ 1978/ pp. 105-113. 

47. Tossell/ J. A.: Calculation of Bond Distances and Cohesive Energies for Gaseous 

Halides Using the Modified Electron Gas Ionic Model. Chem. Phys. Lett., 
vol. 67, no. 2-3/ Nov. 15/ 1979/ pp. 359-364. 

48. Waldman/ Marvin; and Gordon/ Roy G. : Scaled Electron Gas Approximation for 

Intermolecular Forces. J. Chem. Phys./ vol. 71/ no. 3/ Aug. 1979/ 
pp. 1325-1339. 

49. Dalgarno/ A.; DeggeS/ T./ and Williams/ D. A.: Dipole Properties of Molecular 

Nitrogen. Proc. Phys. Soc. (London), vol. 92, pt. 2, no. 576, Oct. 1967/ 
pp. 291-295. 

50. Pauling/ Linus: The Chemical Bond. Cornell Univ. Press, c.1967. 

51. Wigner/ E. P.; and SeitZ/ F. : On the Constitution of Metallic Sodium. II. 

Phys. Rev./ second ser. / vol. 46/ no. 6/ Sept. 15/ 1934/ pp. 509-524. 

52. Brown/ P. J. : A Study of Charge Density in Beryllium. Philos. Mag./ vol. 26, 

no. 6, Dec. 1972, pp. 1377-1394. 

53. Kamaratos, E.; Chang, C. K. ; Wilson, J. W.; and Xu, Y. J. : Valence Bond Effects 

on Mean Excitation Energies for Stopping Power in Metals. Phys. Lett., 
vol. 92A, no. 7/ Nov. 29, 1982, pp. 363-365. 

54. Andersen, H. H.; and Ziegler, J. F.: Hydrogen Stopping Powers and Ranges in All 

Elements. Pergamon Press, Inc., c.1977. 

55. Bakker, C. J.; and Segre, E. : Stopping Power and Energy Loss for Ion Pair Pro- 

duction for 340-Mev Protons. Phys. Rev., second ser., vol. 81, no. 4, Feb. 15, 
1951, pp. 489-492. 

56. Simmons, D. H. : The Range-Energy Relation for Protons in Aluminum. Proc. Phys. 

Soc. (London), vol. 65, pt. A, 1952, pp. 454-456. 

57. Mather, R.; and Segre, E.: Range Energy Relation for 340-Mev Protons. Phys. 

Rev./ second ser., vol. 84, no. 2, Oct. 15, 1951, pp. 191-193. 

58. Sachs, Donald C.; and Richardson, J. Reginald: Mean Excitation Potentials. 

Phys. Rev., second ser., vol. 89, no. 6, Mar. 15, 1953, pp. 1163-1164. 

59. Bogaardt, M.; and Koudijs, B. ; On the Average Excitation Potentials and the 

Range-Energy Relations in Light Elements. Physica, vol. 18, Apr. 1952, 
pp. 249-264. 

60. Vasilevskii, I. M.; and Prokoshkin, Yu. D. : lonizatsionnye Potentsialy Atomov 

(Ionization Potentials of Atoms). Dubna: Ob *edinen. Inst. Yadernykh 
Issledovanni/ 1 960. 


22 



61 


• Shiles, E«; Sasaki, Taizo; Inokuti, Mitio; and Smith, D, Y.: Self-Consistency 

and Sum-Rule Tests in the Kramers-Kronig Analysis of Optical Data: Applica- 

tions to Aluminum. Phys. Rev. B, vol. 22, no. 4, Aug. 15, 1980, pp. 1612-1628. 

62. Andersen, H. H.; Bak, J. F.; Knudsen, H.; and Nielsen, B. R. ; Stopping Power of 

Al, Cu, Ag, and Au for MeV Hydrogen, Helium, and Lithium Ions and Pro- 

portional Deviations From the Bethe Formula. Phys. Rev. A, vol. 16, no. 5, 

Nov. 1977, pp. 1929-1940. 

63. Khandelwal, G. S.: Stopping Power of K and L Electrons. Phys. Rev. A, vol. 26, 

no. 5, Nov. 1982, pp. 2983-2986. 

64. Kahn, David: The Energy Loss of Protons in Metallic Foils and Mica. Phys. Rev., 

second ser., vol. 90, no. 4, May 15, 1953, pp. 503-509. 

65. Neilsen, L. P.: Energy Loss and Straggling of Protons and Deuterons. Mat.-Fys. 

Medd. - K. Dan. Vidensk. Selsk., vol. 33, no. 6, 1961. 

66. Leminens, E.; Fontell, A.; and Bister, M.: Stopping Power of Al, Zn and In for 

0.6-2 .5 MeV Protons. Ann. Acad. Sci. Penn., no. 287, 1968. 

67. Nakata, Hiroshi: Analysis of Energy-Loss Data for 0.2-0. 5 MeV/amu p, a, and N in 

Se. Phys. Rev. B, vol. 3, no. 9, May 1, 1971, pp. 2847-2851. 

68. Sorensen, H.; and Andersen, H. H. : Stopping Power of Al, Cu, Ag, Au, Pb, and U 

for 5-18-MeV Protons and Deuterons. Phys. Rev. B, vol. 8, no. 5, Sept. 1, 

1973, pp. 1854-1863. 


23 


SYMBOLS 


C 


c 

B(t) 

E 

l:(t) 

^Ir^corr 

e 

F(03) 

F (0)) 

D 

f (o)),f 

n 

fn(co) 

f (w) 

P 

g(T) 

H 

h 

tl 

I 


^at 



k 

LPM 


Bohr *s radius 
shell correction 
velocity of light in vacuum 
electric displacement vector 
system excitation energy 
electric field vector 

eigenvalue of nth state 
long-range correlation energy 
magnitude of electron charge 

accumulated oscillator strength with angular frequency below u) 
Dalgarno sum 

optical oscillator strength 

optical oscillator strength of hydrogen atom 

plasma oscillator strength of hydrogen atom 

polarization response 

Hamiltonian operator 

Planck * s cons tant 

h/27T 

mean excitation energy 
atomic mean excitation energy 

molecular mean excitation energy for Bragg's rule with neutral atoms 
mean excitation energy of covalent bonded molecule 

molecular mean excitation energy for Bragg's rule with partially ionic 
atoms 

mean excitation energy of system j 
wave number of ejected electron 
local plasma model 
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L(m) 

I 


m 

n 


R 


AB 


r,ri 


S 

S(m) 

t 


V 


Z 2 ,Z. 

P 


Y 

A 

6 


6 . 

ij 

e(o)) 


\ 

s 

P» Pj^(r ) 
Pn(r) 
p„ (?) 


X 

‘*'0 

CO, (0 


mth moment of log mean excitation energy 

individual electron orbital angular momentum quantum number 

electron mass 

density of atoms 

stochiometric coefficient 

Rydberg's constant 

internuclear distance in molecule composed of atoms A and B 
electron position vector 
radius of Wigner-Seitz sphere 
stopping power 

mth moment of oscillator strengths 
time 

ion velocity 
ion charge number 
atomic number 

ratio of ion velocity to velocity of light 
plasma frequency shift 
straggling parameter 
polarization correction 

chemical bond shift in mean excitation energy 
dielectric function 
average electron separation 
electron density 

electron density in hydrogen atom 
electron density in hydrogen molecule 
screening parameter 

ground state electronic wave function 
angular frequency of oscillator 


classical plasma frequency 


0) 

P 

< > averaged value 

An arrow over a symbol denotes a vector, A prime over a symbol denotes a 
variable of integration. 


26 


































TABLE III.- MOLECULAR MEAN EXCITATION ENERGY 


Chemi cal 
species 

Molecular mean excitation energy, 

eV, for - 

Present 

theory 

Experimental 

empirical 

Bragg* s 
rule 

CH4 

44.7 

42.8 

35.1 

<CH2>x 

55.0 

53.4 

43.5 

C6»6 

60.6 

61.4 + 1.9 

50.6 

«2 

18.9 

18.5 ± 0.5 

15.0 

Graphite 

76.1 

78.5 ± 1.5 

62.0 


TABLE IV.- MOLECULAR MEAN EXCITATION ENERGIES FOR COVALENT GASES 


Chemical 

species 

R , a 
AB 0 

I, eV 
(a ) 

I, eV 
(b) 

»2 

1 .40 

1 8.9 

19.2 ± 0.4^ 


2.08 

85.0 

82.0 ± 1 . 6 ^ 

O2 

2.34 

99.6 

95 ± 1.9^^ 

^2 

2.67 

114.2 

115 + 10 

C ^2 

3.76 

170.8 

171 ±14 


Local plasma model. 

^Reference 40. 

These values are strongly influenced by reference 11 



















TABLE V.- MEAN EXCITATION ENERGIES FOR COVALENT-BONDED CRYSTALS 


Chemical 

species 

R , a 
AB o 

I, eV 
(a) 

I, eV 
(b) 

B (tetragonal) 

3,06 

67.3 

76 ± 7.6 

C (diamond) 

2,94 

75.3 


C (graphite) 

2,68 

76.1 

78 ± 2.3 

Si (diamond) 

4.42 

151.0 

173 ± 4 

P (black) 

4.16 

155.7 

181 ± 14^ 

S (rhombic) 

3.85 

162.7 

190 ± 15^ 


Local plasma model. 
Reference 40, 

Unspecified allotropic form. 


TABLE VI,- IONIC-BOND PARAMETERS 


Chemical 

species 


P 

M 

(0 

< 

IlB, eV 

Ig, eV 

I, e 

0.50 

97.6 

91.7 

91,0 

96.4 

.25 

27.8 

25.2 

25.9 

26.7 

.90 

83.4 

92.6 

81.6 

93.6 
























TABLE VII 


MEAN EXCITATION ENERGIES OF IONIC CRYSTALS 


Chemical 

species 

^AB' 

P 

I, eV 
(a) 

I, eV 
(b) 


3.80 

0.90 

92.8 

94 

± 8 

LiCJl 

4.86 

.73 

139.1 

144 

± 12 

NaF 

4.37 

.91^^ 

131.5 

147 

± 12 

NaCJl 

5.31 

.75 

159.1 

181 

± 14 


Local plasma model, 

^Reference 40. 

^Pauling partial ionic character function. 


TABLE VIII.- METALLIC PARAMETERS FOR SELECTED METALS OF FIRST TWO ROWS 


Chemical 

species 

(a) 

s o 

Lithium 

34.0 

3.260 

Be ry Ilium 

38.6 

2.375 

Sodium 

123.6 

3.99 

Magnesium 

121.2 

3.34 

Aluminum 

124.3 

2.991 


‘=^Ref erence 6. 
^^igner-Seitz model 
^Reference 40. 


41 .5 ± 3.7 
63.7 ± 3.2 
162 ± 8 
164 ± 8 
166 ± 3 
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Figure 1 •- Cumulative oscillator strength distribution 
for atomic hydrogen and helium. 
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Figure 2.- Moments of oscillator strengths of hydrogen atom 
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Figure 5.- Comparison of atomic mean excitation 
energies according to quantum calculations of 
Inokuti et al. and local plasma model. Empirical 
values are from reference 40. 
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Figure 6.- Comparison of atomic straggling parameters 
according to quantum calculations of Inokuti et al 
(refs. 6 to 8) and local plasma model. 






Figure 7.- Moments of oscillator strengths from empirical 
values of reference 49 compared with local plasma model 
using Gordon-Kim molecular model densities. 
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Figure 8.- Mean excitation energies for 
partially ionic atoms. 
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Figure 9.- Comparisons of reduced stopping power for 
aluminum for several mean excitation energies with 
range of experimental data. (See refs. 60 and 62 
to 66. ) Shaded area is band of experimental data. 
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^Atomic values pi LPt* 
ODahmer et at. (refs. 6 to 8J 
aMolecular values by LPM 
ASeltzer-Berger [ref. 40 J 
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Figure 10,- Mean excitation energies for atoms, molecules, 
solids, and metals. Specific data taken from tables III 
to V and table VIII. 
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